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The vortex states in a thin mesoscopic disk are investigated within the phenomenological 
Ginzburg-Landau theory in the presence of a step-like external magnetic field with zero average 
which could model the field resulting from a ferromagnetic disk or a current carrying loop. The 
regions of existence of the multi-vortex state and the giant vortex state are found. We analyzed 
the phase transitions between them and found regions in which ring-shaped vortices contribute. 
Furthermore, we found a vortex state consisting of a central giant vortex surrounded by a collection 
of anti- vortices that are located in a ring around this giant vortex. 



I. INTRODUCTION 



Recent progress in microfabrication and measurement techniques makes it possible to study the properties of 
mesoscopic superconducting samples, and, in particular, disks with sizes comparable to the penetration depth A and 
the coherence length £ |l]-^| . Depending on the disk size various superconducting states can exist in such samples in 
an external magnetic field. 

In small disks, where the confinement dominates, there exist only giant vortex (GV) states (circular symmetric 
states with a fixed value of angular momentum L) ^|J^] . For sufficiently large disks the GV state can break up into 
the multi- vortex (MV) state (a mixture of giant vortices with different angular momentum) HQ. Also ring-shaped 
vortex states with larger energy than GV and MV ones were predicted Q . These ring-shaped two-dimensional vortex 
states have a cylindrically symmetric magnetic field profile and they are different from the ring-vortices that were e.g. 
found in three-dimensional superfluid liquid helium. But the analysis of Ref. Q showed that in a superconducting 
disk in the presence of a homogeneous magnetic field the ring-shaped vortices are unstable. 

In the present paper we study the MV states of a thin superconducting disk in an inhomogeneous magnetic field 
Hq with step-like profile. This magnetic field profile could approximate the field resulting from a magnetic dot (they 
are usually used in experiments, e.g. see Ref. ||) as well as the field of a current loop (Fig. ^). We found that this 
model field profile captures a lot of interesting physics and that the step height and width strongly influence the 
superconducting phase diagram. 
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FIG. I. The magnetic field profile from a magnetic disk (a) and a current loop (b) and the corresponding step field models. 



II. APPROACH USED 

The details of our approach can be found in Refs. ||,||,^| and, therefore, we briefly sketch only the necessary steps. 
We restrict ourselves to sufficiently thin disks such that d << £, A. In this case, to a first approximation, the magnetic 
field due to the circulating superconducting currents may be neglected and the magnetic field inside the disk equals 



1 



the external one. Following Refs. |3|-|6|], the order parameter of the MV state is written as a linear combination of 
eigenfunctions of the linearized Ginzburg-Landau equation 



exp (iLjcp) , 



(1) 



L i= n=0 



where p is the radial distance from the disk center, ip is the azimuthal angle, and n enumerates the different radial 
states for the same Lj . In the homogeneous magnetic field case L is the value of the effective total angular momentum 
which is equal to the number of vortices in the disk. Later we will see that, for our inhomogeneous magnetic field 
case, the assignment of the total vorticity can be tricky. 

Substituting Eq. ([I]) in the free energy (measured in F = H^V/8tt units) 



F = T 



dV 



> 4 
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(here h (p) = V x A(p) is the local magnetic field, A is the vector potential) we obtain F as a function of the 
complex parameters C n j Ji . Minimization of F with respect to these parameters allows to find the equilibrium vortex 
configurations and to determine their stability. 

We restrict ourselves to states built up by only two components in Eq. (Q). It brings quantitative bounds in our 
analysis but, nevertheless, gives the correct qualitative behavior and facilitates the physical insight into the problem. 
Using two different approaches from Refs. (||| AfV configurations with three and five components in Eq. (|l|) were 
checked. This resulted in 1) a reduction of the region of existence of metastable MV states at the low magnetic field 
limit, 2) an appearance of additional unstable states corresponding to the saddle points of the free energy function, 
and 3) no new vortex configurations in the ground state. Thus, two components in Eq. (|lj) are sufficient to describe the 
vortex structure in the ground state of our system. Depending on the radial state of the components we defined two 
different configurations: giant-giant (m = ri2 = 0) and giant-ring [n\ = 0, ri2 = 1) MV states. The first represents a 
combination of two GV states while the second is composed of one giant and one ring vortex. 



III. CENTERED STEP-LIKE MAGNETIC FIELD 



First, we investigate the case of a centered step-like magnetic field corresponding to a magnetic disk. The effect 
of the size of the superconducting disk on the vortex state is shown in the phase diagram of Fig. ^|. The parameters 
considered are R\ = 4.5£ and i?2 = 6.0£. The solid lines indicate where the ground state of the free energy changes 
from one state to another and dashed lines correspond to transitions between different GV states as metastable states. 
The right axis of the diagram shows the corresponding positive flux through the region p < Ri . 
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FIG. 2. The Hi n — R phase diagram for the ground state of a thin superconducting disk in a step-like magnetic field. Solid 
curves indicate transitions between different vortex states including different MV regions. Dashed curves denote the transitions 
between different metastable GV states. 
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One can clearly see the re-entrant behavior. For example, for R = 9.0£ we observe the following change of total 
vorticity L — — >1^2^3^1^0. Notice that the ground state with total vorticity L = and L = 1 covers the 
largest part of the phase diagram. With increasing disk size, all other GV states are strongly suppressed in favor of 
the various MV states. For small disk sizes different GV states are present as the ground state. But with increase of 
disk size, islands with different MV configurations (giant-ring and giant-giant MV states) dominate the ground state 
diagram. For R — 6.38£, a giant-ring MV state appears, and for R = 6.7£ we obtain giant-giant states. 

As one can see in Fig. ||, for example, for R = 9.0£ and with increasing magnetic field Hi n , giant-giant combinations 
(Li : L 2 ) = (0 : 7), (0 : 8), (0 : 9) etc. become the ground state. In the contour plot of the | ^-distribution for these 
states we observe single vortices which are arranged on a ring with a low density area (see Fig. ||(a), blue color - zero 
density, red - highest density) situated in the center of the disk (see Fig. |](a), for the (0 : 7) state). This central area 
is associated with a giant vortex, and encircling it leads to a phase change showing the vorticity L = 7. 
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FIG. 3. (a) Contour plot of the superconducting density for the (0 : 7) giant-multi anti-vortex ground state, and (b) the 
corresponding contour plot of the phase. Notice that the phase near the boundary is near (blue) or 2n (red) but due to the 
finite numerical accuracy it oscillates between 2n — e and 2ir + e where e « 10 -5 . 

Fig. [|(b) shows the contour plot of the corresponding phase with seven anti-vortices arranged in a circle around 
the giant vortex and the total vorticity is L = 0. It is important to emphasize that for the disk size R = i?2 and 
with changing of R\ giant-giant MV states consist only of vortices arranged radially symmetric with no anti-vortices 
present. Thus, it is obvious that enlarging the superconducting disk enhances the influence of the negative part of 
the step-like magnetic field. As a result we obtain a new giant-giant form called giant-multi anti-vortex states. The 
total vorticity is now equal to the lowest vorticity of the giant vortex states of which the multi-vortex consists of. 
This inevitably leads to a re-entrant behavior in the total vorticity, which is in agreement with Ref. |]lof where it was 
found numerically that several vortices can enter (or exit) at once for disks with sufficiently large radius. For R = R^ 
disks we only found that L — ► L ± 1 transitions are possible. With increase of the disk radius the (0 : 8) and (0 : 9) 
giant-multi anti-vortex states become (1 : 8) and (1 : 9), respectively. For R > 12.33£, the multi-vortex states become 
metastable again and the L = and L = 1 states are now the only ground states. Moreover, for R > 18.67£ we have 
the Meissner state for all values of the applied magnetic field. 

IV. SHIFTED STEP-LIKE MAGNETIC FIELD 

Next, we move the step-like field profile along the radius of the superconducting disk, i.e. consider a ring magnetic 
field profile as in the case of a current loop. In Fig. ||, we present the phase diagram as function of Rd- The parameters 
of the magnetic field profile are R\ ~ Rd = 4.5£, R2 — Rd = 6.0£, and R = 9.0£. Thick solid curves indicate the 
transitions between different vortex states and dashed lines denote transitions between metastable GV states. 
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FIG. 4. The Hi n — Rd phase diagram for the ground state of a thin superconducting disk in a ring step-like magnetic field 
distribution. The same curve convention is used as in Fig. ^. 

The re-entrant behavior, as in previous case, is clearly visible. The most important result is that the MV states 
dominate this diagram. Moreover, we have both types of multi-vortices, i.e. giant-giant and giant-ring, as ground 
state, and with shifting the field profile towards the disk periphery, giant-giant MV configurations cover most of the 
superconducting region. The giant-giant MV states in this case appear both as states with no anti-vortices present, 
and as giant-multi anti-vortex states. Following the transition lines, a correspondence between different giant-multi 
anti- vortex states {L\ : L2) can be seen: L2 remains the same, while L\ increases from to 2. As one can see, the 
latter is strongly correlated with re-entrant behavior. In the rest of the diagram, the MV states with the "classical" 
geometry (both giant-giant and giant-ring) dominate. However, a difference between those states exists. As shown in 
Fig. |^, in the case of a giant-ring state, vortices are preferentially distributed within a ring-shaped lower density area. 
Further, considering the high density areas at the disk periphery, a shift in phase of A8 = tt/L is observed, where L 
is the total vorticity. Although these states have a different origin sometimes they can exhibit a similar distribution 
of vortices but the phase is always able to distinguish between them. 
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FIG. 5. Contour plots of the superconducting wave density for giant-giant (a) and giant-ring (b) MV state with the total 
vorticity L = 4. 
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